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We present a theoretical study of the ground state of the BCS-BEC crossover in dilute two-dimensional Fermi 
gases. While the mean-field theory provides a simple and analytical equation of state, the pressure is equal to 
that of a noninteracting Fermi gas in the entire BCS-BEC crossover, which is not consistent with the features of 
a weakly interacting Bose condensate in the BEC limit and a weakly interacting Fermi liquid in the BCS limit. 

The inadequacy of the 2D mean-field theory indicates that the quantum fluctuations are much more pronounced 
than those in 3D. In this work, we show that the inclusion of the Gaussian quantum fluctuations naturally 
recovers the above features in both the BEC and the BCS limits. In the BEC limit, the missing logarithmic 
dependence on the boson chemical potential is recovered by the quantum fluctuations. Near the quantum phase 
transition from the vacuum to the BEC phase, we compare our equation of state with the known grand canonical 
equation of state of 2D Bose gases and determine the ratio of the composite boson scattering length ag to the 
fermion scattering length a 2 D- We find ag - 0.56fl2D, in good agreement with the exact four-body calculation. 

We compare our equation of state in the BCS-BEC crossover with recent results from the quantum Monte Carlo 
simulations and the experimental measurements and find good agreements. 

PACS numbers: 05.30.Fk, 03.75.Ss, 67.85.Lm, 74.20.Fg 


I. INTRODUCTION 

The experimental realization of ultracold atomic Fermi 
gases with tunable interatomic interactions has opened a new 
era for the study of some longstanding theoretical proposals 
in many-fermion systems. One interesting proposal is the 
smooth crossover from a Bardeen-Cooper-Schrieffer (BCS) 
superfluid ground state with largely overlapping Cooper pairs 
to a Bose-Einstein condensate (BEC) of tightly bound bosonic 
molecules - a phenomenon suggested many years ago M- 
A simple but important system is a dilute attractive Eermi gas 
in three dimensions (3D), where the effective range of the 
short-ranged interaction is much smaller than the interparti¬ 
cle distance. The system can be characterized by a dimen¬ 
sionless gas parameter l/ikya^u), where aso is the i-wave 
scattering length of the short-ranged interaction and kp is the 
Eermi momentum in the absence of interaction. The BCS- 
BEC crossover occurs when the parameter 1 / (kpa^u) is tuned 
from negative to positive values Hi], and the BCS and 
BEC limits correspond to the cases lUkpaio) and 

ll(kpa^Y)) —> respectively. 

The BCS-BEC crossover phenomenon in 3D dilute Eermi 
gases has been experimentally demonstrated by using ultra¬ 
cold gases of ®Li and ""'K atoms iHHH, where the i-wave 
scattering length and hence the gas parameter l/ikpa^n) were 
tuned by means of the Eeshbach resonance lll2.fo . The equa¬ 
tion of state and various static and dynamic properties of the 
BCS-BEC crossover have become a big challenge for quan¬ 
tum many-body theory ifldl - l^ because the conventional per¬ 
turbation theory is no longer valid. At the so-called unitary 
point where > °o, the only length scale of the system 

is the inter-particle distance. Therefore, the properties of the 
system at the unitary point IKkpa^u) - 0 become universal, 
i.e., independent of the details of the interactions. All thermo¬ 


dynamic quantities, scaled by their counterparts for the non¬ 
interacting Eermi gases, become universal constants. Deter¬ 
mining these universal constants has been one of the most in¬ 
triguing topics in the research of the cold Eermi gases ll25l - [^ . 

On the other hand, it was suggested that a 2D Eermi gas 
with short-ranged i-wave attraction can also undergo a BCS- 
BEC crossover Unlike 3D, a two-body bound state 

always exists in 2D even though the attraction is arbitrarily 
weak. The BCS-BEC crossover in 2D can be realized by 
tuning the binding energy of the bound state. Studying the 
BCS-BEC crossover in 2D will help us understand the physics 
of pseudogap and Berezinskii-Kosterlitz-Thouless transitions 
in fermionic systems 1^ . In recent years, quasi-2D atomic 
Eermi gases have been experimentally realized and studied 
by a number of groups Bbl - mn . In cold-atom experiments, 
a quasi-2D Eermi gas can be realized by arranging a one¬ 
dimensional optical lattice along the axial direction and a 
weak harmonic trapping potential in the radial plane, such 
that fermions are strongly confined along the axial direction 
and form a series of pancake-shaped quasi-2D clouds. The 
strong anisotropy of the trapping potentials, namely, Wj » a>± 
where co^ (wj.) is the axial (radial) frequency, allows us to use 
an effective 2D Hamiltonian to deal with the radial degrees of 
freedom. Experimental studies of quasi-2D Eermi gases have 
promoted great theoretical interests in the past few years lld^ 

It is known that in 3D, even the mean-field theory predicts 
that the system is a weakly interacting Bose condensate in 
the strong attraction limit The composite boson scatter¬ 
ing length is shown to be ab = The inclusion of 

Gaussian pair fluctuations iflbl - fisl] recovers the Eermi liquid 
corrections in the weak attraction limit and modifies the com¬ 
posite boson scattering length to^^ 0.55a3D, which is close 
to the exact result - O.baan Moreover, the equation 
of state (EOS) in the BCS-BEC crossover agrees excellently 
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with the quantum Monte Carlo results and the experimental 
measurements if the Gaussian pair fluctuations are taken into 
account CMl. In contrast, the mean-field theory for 2D 
Fermi gases does not predict a weakly interacting 2D Bose 
condensate in the strong attraction limit ll^ 1^ . The cou¬ 
pling constant between the composite bosons is predicted to 
be energy independent, which arises from the inadequacy of 
the Born approximation for four-body scattering in 2D. As a 
result, the 2D mean-field theory predicts that the pressure of a 
homogeneous 2D Fermi gas is equal to that of a noninteract¬ 
ing Fermi gas in the entire BCS-BEC crossover. However, re¬ 
cent experimental measurements llj^lddT] and quantum Monte 
Carlo simulations ll4bl - l4^ show that the pressure in the strong 
attraction limit is vanishingly small in comparison to that of a 
noninteracting Fermi gas, which is consistent with the picture 
that the system is a weakly interacting 2D Bose condensate. 
These results indicate that the 2D mean-field theory is not ad¬ 
equate even at the qualitative level, and quantum fluctuations 
are much more important in 2D. 

In analogy to the 3D case, we expect that the inclusion of 
Gaussian pair fluctuations in 2D naturally recovers the feature 
of the weakly interacting 2D Bose condensate in the strong 
attraction limit. This has been demonstrated recently by using 
the pole approximation for the Goldstone mode and the di¬ 
mensional regularization for the untraviolet divergence, which 
leads to an eleg ant derivation of the composite boson scatter¬ 
ing length 117^ . However, the pole approximation is limited 
in the strong attraction limit because of the use of the Bogoli- 
ubov dispersion for the Goldstone mode. We also note that 
the ultraviolet (UV) divergence arising from the pole approx¬ 
imation can be naturally avoided in the full treatment of the 
collective modes fla - USl] . In this work, we study the influ¬ 
ence of quantum fluctuations on the EOS of 2D Eermi gases in 
the entire BCS-BEC crossover. With the full EOS beyond the 
pole approximation, we determine the ratio of the composite 
boson scattering length to the fermion scattering length 020 
by comparing our EOS near the vacuum-BEC quantum phase 
transition with the known grand canonical EOS of weakly in¬ 
teracting 2D Bose gases We obtain as - 0.56a2D, in 

good agreement with the exact four-body calculation and 
the pole approximation treatment if?^ . We also perform nu¬ 
merical calculations for the canonical EOS of a homogeneous 
2D Eermi gas in the BCS-BEC crossover. In addition to re¬ 
covering the weakly interacting Bose condensate in the st rong 
attraction limit, we find that the Eermi liquid corrections 11831 - 
[sH] can also be recovered at sufficiently weak attraction. We 
compare our EOS with the recent results from quantum Monte 
Carlo simulations and experimental measurements in the en¬ 
tire BCS-BEC crossover and find good agreements. 

The paper is organized as follows. We set up our theoretical 
framework for 2D Eermi gases beyond mean field in Sec. HI] 
We study the strong attraction limit and determine the com¬ 
posite boson scattering length in Sec. |III| We present our the¬ 
oretical predictions for the EOS in the BCS-BEC crossover 
and compare our results with the quantum Monte Carlo data 
and experimental measurements in Sec. |IV] We summarize in 
Sec. |V] The natural units ^ = 1 are used throughout. 


II. HAMILTONIAN AND GRAND POTENTIAL 

We consider a spin-1/2 (two-component) Eermi gas in two 
spatial dimensions with a short-ranged s-wave attractive inter¬ 
action between the unlike spins. In the dilute limit, the inter¬ 
action potential can be safely modeled by a contact interac¬ 
tion. The grand canonical Hamiltonian density of the system 
is given by 

H= ^ i/ro-(r)'7/oiAtr(r)-t/iAT(r)iAi(r)'Ai(r)'AT(r)> (1) 

o'=T,i 

where i/''f(r) and lAj(r) represent the annihilation field opera¬ 
tors for the two spin states of fermions, 'Ko = -V^/(2m) -p is 
the free single-particle Hamiltonian, with m being the fermion 
mass and fi being the chemical potential, and U > 0 de¬ 
notes the i-wave attractive interaction occurring between un¬ 
like spins. 

The contact coupling U is convenient for performing the¬ 
oretical derivations. However, it should be renormalized by 
using some physical quantities so that we can obtain finite re¬ 
sults in the many-body calculations. With the contact interac¬ 
tion U, the Lippmann-Schwinger equation for the two-body T 
matrix reads 

Tll{E)^-U-^-n{E), (2) 

where E - I mis the scattering energy in the center-of-mass 
frame and the two-particle bubble function n(£’) is given by 

n(E) = y , ■ (3) 

^ E + ie- 2ek 

Here e = 0^ and Sk = k^/(2m). We use the notation 
2k = / £/“k/(27r)^ throughout. The cost of the use of the 
contact interaction is that the integral over k suffers from UV 
divergence. We regularize the UV divergence by introducing 
a hard cutoff A for |k|. Eor large A we obtain 

m A^ 

n(£) =-In — -H — In (-£ - ie) . (4) 

4;r m 47r 

Next we match the scattering amplitude f{k) - {An jm)T 2 b{E) 
to the known 2D s-wave scattering amplitude in the zero- 
range limit, which is given by f{k) — l/[ln(£B/E) + i^] s 
113] ■ Here sb is the binding energy of the two-body bound state 
which characterizes the attraction strength. Thus we obtain 

1 ^ In = y 1 (5) 

t/(A) 47r ” msB |j^ 2sk + sb ' 

The above results should be understood in the large-A limit. 
After the renormalization of the bare coupling U through the 
physical binding energy sb, the UV divergence in the many- 
body calculations can be eliminated and we can set A —> 00 to 
obtain the final finite results. 

In the imaginary-time functional path integral formalism, 
the partition function at temperature T is 

r[u!(A][#]exp{-.S[^,^]}, (6) 
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where the action 


A. Mean-fleld approximation 


*S[iA, *A] = J dx [il/driff + Hit//, if/)]. (7) 

Here x = (r, r), with r being the imaginary time, and J dx - 

dr J d~r with /3 - l/T. To decouple the interaction term, 
we introduce an auxiliary pairing field <l)(;r) and apply the 
Hubbard-Stratonovich transformation. Then the the partition 
function can be expressed as 


Z = exp <[),<[)*]}, (8) 


where the action reads 

Q I ^ 

o = I ax 


= J" dx -—- J" dx J" dx^{x)G ^{x,x)^{x). (9) 


Here the Nambu-Gor’kov spinor i/^(x) = {^'\{x),^i{x)Y is 
employed and the inverse Green’s function G^^{x,x') in the 
Nambu-Gor’kov representation is given by 


G *(x, x') = 


-dr - 'Ho 

(fi*(x) 


<l>(x) 

-dr+^0 


jb(x - x')- 


( 10 ) 


Integrating out the fermion fields, we obtain 

Z = J[ufO][of(I)*]exp{-^eff[0,0*]}, (11) 

where the effective action reads 

>SefF[0, J- Trln[G-'(x, x')]. (12) 

The partition function cannot be evaluated analytically 
since the path integral over <t) and €>* cannot be carried out. 
At r = 0, the pairing field 0(x) acquires a nonzero and uni¬ 
form expectation value (<l>(x)) = A, which serves as the or¬ 
der parameter of superfluidity. Due to the U(l) symmetry, we 
can set A to be real and positive without loss of generality. 
Then we write <fi(x) = A H- <pix), where <p{x) is the fluctuation 
around the mean field. The effective action >Sefr[<b, <b*] can be 
expanded in powers of the fluctuation 0(x); that is. 


>Seff[0, O*] = *Smf + 01 + ■ ■ ■ , (13) 

where >Smf = ‘5efF[A, A] is the saddle-point or mean-field 
effective action and >Sgf[ 0, 0*] is the contribution from the 
Gaussian fluctuations (GFs). The higher-order contributions 
from non-Gaussian fluctuations are not shown. Accordingly, 
the grand potential can be expressed as 

G = -F Ogf + ■ ■ ■ » (14) 


In the mean-field approximation, the effective action is ap¬ 
proximated as >SefF[<5, <5*] - 'Smf- The quantum fluctuations 
are completely neglected. At T =0, the mean-field grand 
potential can be evaluated as 

A2 

Gmf = — + ^ (^k - T^k) , (15) 


where = Sk - and isk = + A^. The UV divergence 

can be eliminated by using Eq. ©. We obtain 

QMF = A^Vf-— - V (16) 

■'^\2Sk + SB Ek+^k/ 

The order parameter A should be determined as a function of 
ju by using the extreme condition (9 Qmf/ 0A = 0. We obtain 
the gap equation 


1 

U 



(17) 


or, explicitly. 



1 

2ek + SB 



(18) 


It is very fortunate that in 2D the integral over k can be 
carried out. The grand potential reads 


Gmf — 


mA^ 

4-k 


In 




1 


.(19) 


Using the extreme condition (9 Qmf/ 0A = 0, we obtain 

Amf)/^) = VeB(2/t H- en) ©(2// -F sb), (20) 


which determines analytically the order parameter A as a func¬ 
tion of the chemical potential yU. Substituting this result into 
Qmf, we obtain the mean-field grand canonical EOS 

Gmf(/2) = -^(2jU -F Ssfeil/i + sb). (21) 

on 

The mean-field contribution to the particle density is given by 
m 

«MF(/t) = :^('2n + eB)©(2jU -f sb)- (22) 

2n 

The above mean-field results show that the system undergoes 
a second-order quantum phase transition from the vacuum to a 
matter phase with nonvanishing density. The critical chemical 
potential is given by 



where Gmf = ‘5mf/(/ 1V), with V being the volume, and Ggf 
is the contribution from the Gaussian fluctuations. 


( 23 ) 
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B. Gaussian pair fluctuation theory 


Now let us include the quantum fluctuations. We include 
the Gaussian fluctuations only and approximate the effective 
action as >Seff[0,O*] - *Smf + '5gf[ 0,0*]- The advantage 
of this Gaussian approximation is that the path integral over 
0 and (ff can be carried out analytically. To evaluate the 
quadratic term >Sgf[ 0, <^*], we make the Fourier transforma¬ 
tion 

m ^ (24) 

Q 

where Q - {iqi, q), with qi = llnT (/ 6 Z) being the boson 
Matsubara frequency. We use the notation YjQ - TYjiTjq 
throughout. After some manipulations, >Sgf[<^, 0*] can be ex¬ 
pressed in a compact form, 

>Sgf[<^, ) M(0 J ) (25) 

The inverse boson propagator M(2) is a 2 x 2 matrix. At 
r = 0 , its elements are analytically given by 


Mn( 0 =M 22 (-e) 


1 + y 

u 4-‘ 




2 2 


^ / £k ^k+q ^ ^/ ■*" ^k "I" ^k+q 


= z 


Mi2(0=M2i(0 

^k^k^^k+q^k+q 


Mk^'k^k+q^k+q 
iqi -H Fk + Fk+q iqi - £k - £k+q 


(26) 


Here the BCS distribution functions are defined as i/^ = (1 - 
fk/Fk)/2 and u^-1- u^. Note that Eq. (|3 should be used to 
eliminate the UV divergence. 

Considering the Gaussian fluctuations only, the partition 
function is approximated as 

Z - exp(-*SMF) Jexp{ - .Sgf[0,<^*]}. (27) 

Carrying out the path integral over (ft and (jf, we obtain the 
grand potential G = Gmf + ^^gf, where the contribution from 
the Gaussian fluctuations can be formally expressed as 


where the phase shifts are defined as bM(w, q) = 
-Im In det M(w + /e, q), (5ii(w, q) = -ImlnMii(w -i- fe, q), 
and 622 {to, q) = -ImlnM 22 (w -i- ie, q). 

A crucial element of the Gaussian pair fluctuation (GPF) 
theory is that the order parameter A should be determined by 
the extreme of the mean-field grand potential f^F rather than 
the full grand potential Q. = Qmf + Ggf ifl^fT^ . Therefore, 
we still use the mean-field gap equation or the analytical re¬ 
sult, (I 2 TI 1 . The advantages of the use of the mean-field gap 
equation can be summarized as follows: 

(i) The mean-field solution for the order parameter, (|2T]) . guar¬ 
antees Goldstone’s theorem. The dispersion of the Goldstone 
mode can be obtained by solving the equation 

detM(w,q) = 0 (31) 

for a> smaller than the two-particle continuum. The use of the 
mean-field solution, (l2T]) . for the order parameter ensures that 
detM(0,0) = 0. Therefore, the lightest collective mode is 
gapless and has a linear dispersion at low momentum q. We 
expect that the most important contribution from the quantum 
fluctuations is the Goldstone mode fluctuation. The use of the 
mean-field gap equation ensures that the Goldstone mode is 
gapless and hence enables us to take into account correctly 
the contribution from the Goldstone mode. 

(ii) The use of the mean-field solution, (|2T]) . maintains the fa¬ 

mous Silver Blaze property lls^ Is^ even if we consider the 
contributions from the quantum fluctuations. Even though the 
critical chemical potential jUc = -£ 3/2 for the vacuum-matter 
transition is obtained from the mean-field approximation, we 
expect that it is exact because the minimal chemical potential 
to create a bound state is exactly 2 / 2 ^ = -en- Eor fi < jj.c and 
at r = 0 , the system stays in the vacuum phase with vanishing 
pressure and density. This is known as the Silver Blaze prob¬ 
lem Obviously, the mean-field EOS satisfies this 

property. Now we show that the Gaussian contribution Ogf 
also satisfies this property. Eor fj. < fj.c, we have A = 0 and 
hence Ogf is given by 

QGF = yinMo(/^,,q)e'4'‘’", (32) 

Q 

where the pair susceptibility Mo(f^/, q) in the vacuum phase 
is analytically given by 


Ogf = lndetM(0. 

G 


(28) 


However, this formal expression is divergent because the con¬ 
vergent factors are not appropriately considered. Considering 
the convergent factors leads to a finite result EQ: 



Mii(0 

M 22 (e) 


detM(0 




(29) 


The Matsubara frequency sum can be converted to a standard 
contour integral. At T = 0, we have 



diL) 


[bM(w, q) -F bii(w, q) - b 22 (w, q)] ,(30) 


Mo(/^;,q) = ^(7 


1 


1 


i?/ “ fk+q/2 - ^k-q/2 2ek + 


m 

= 

4n 


SB 


At r = 0 , we obtain 
f^GF = 




bo(w, q), 


(33) 


(34) 


where 60 (a), q) = -ImlnMo(w + ie, q). It is easy to show that 
6o(o), q) = 0 for cu < 0 in the vacuum phase ^ Therefore, 
we have exactly Qgf = 0 for ji < Accordingly, the particle 
density also vanishes in the vacuum. 
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C. Imaginary frequency integration formalism 


For the Gaussian contribution Qgf, it is convenient to em¬ 
ploy an alternative formalism which automatically satishes 
the Silver Blaze property and also leads to faster convergence 
for numerical calculations. To this end, we dehne two func¬ 
tions, M^j(z, q) and q) 1131 , which are given by 




U 


z 




Ek - Ek 


(35) 


+q 


Using the gap equation, (O, and the fact that < 1, we 
can show that q) has no singularities or zeros in the 

left half-plane (Rez < 0). Therefore, the Matsubara sum 


Yiq, lnM^j(i^;, q) vanishes at T = 0 since lnM^j(z, q) has no 
singularities in the left-half plane. Therefore, the Gaussian 
contribution at T = 0 can be expressed as Cl 



Muiiqi, q)M 22 {iqi, q) - M^^iiqi, q) 


■ (36) 


At r = 0, we replace the discrete Matsubara frequency sum 
with a continuous integral over an imaginary frequency; i.e.. 


OO ^CO 1 

Tj^Xiiq,)^ J (37) 

l—~oo 

After some manipulations, we obtain 


J 




In 


1 oa4/ A(oJ,q)C(oj,q) + OJ^B{a),q)D(aj,q) + 2F^(aj,q) g C^(aj,q) + a/D^(uj,q) 
1 - 2A i^i) - -J- ------H- A (//)■ 


1 1 


iEk+ql2 + ^k+q/2)(£’k-q/2 + ^k-q/ 2 ) 

2ek + eB 4 

,£’k+q/2 Ek-ql2i 

(£’k+q/2 + £’k-q/2)^ + 


A^{a), q) -H q) A^{co, q) -h to^B^ito, q) 

where we have used the fact that the integrand is real and even in o). The functions A, B, C, D, and F are dehned as 

A(w,q) = Yj 

k 

B(w,q) = 2 

k 

C(cj,q)^Y 

k 

D{u,q)^Y 


1 


{Ek+q/2 + ^k+q/2)(£’k-q/2 + ^k-q/ 2 ) 


4£'k+q/2£’k-q/2 (£’k+q/2 + Ek-q/ 2 )^ + Ci? 

1/1 1 \ 1 

-I- 


1 


4 (Tik+q/Z £’k-q/2 / (£’k+q/2 + ^k+q/2)(£’k-q/2 + ^k-q/2) (Uk+q/2 + £’k-q/2)^ + 
1 1 


4£’k+q/2£’k-q/2(£’k+q/2 + ^k+q/2)(£’k-q/2 + ^k-q/ 2 ) (Ek+q/l + Ek-q/l)^ + O)^ ’ 


U(u2,q) = ^- 


1 / 1 


1 


1 


4 \£'k+q/2 Ek~ql2j (£’k+q/2 + £’k-q/2)^ + 

I- 


(38) 


(39) 


Note that A(/2) is given by the mean-held solution, (|2TI) . The 
BCS-type dispersions £’k±q/2 in ( l39] l are hence analytically 
given by 

£'k±q/2 = -y/(ek±q/2 “ + Sb(2zU -H SB)0(2yU -H Sb)- (40) 

The integrand in (l38T l vanishes in the vacuum jj. < fic and 
hence the Silver Blaze property is automatically satished. 
Moreover, because we use the mean-held gap equation, (|2TI) . 
we can replace 1 / {2sk -i- sb) with 1 / (2Ek) in the expression of 
the function A(w, q). Then we hnd that the integrand in (l38l l 
diverges near (o), q) = (0,0), which indicates that the most im¬ 
portant contribution is from the low-energy Goldstone mode. 

In summary, the grand canonical EOS in the GPF theory is 
given by 

Q(fj) - Qmf(a) + (41) 

The particle density «(//) reads 

n(fj) = riuFifj) + ncFifE), (42) 


where the GF contribution is formally given by 

d^lcFiM) 


ncFifj) 


dfx 


Here dldji represents the full derivative with respect to /r; 


dElcFilJ) 

dfi 


^Qgp ^Qgp dA 

-F 


dfi 


dA dfi 


(43) 
i.e., 

(44) 


In 3D, it was shown that the second term is crucial to produce 
in the BEC limit the composite boson scattering length - 
0.55a3D which is very close to the result ab = 0.6^30 

from the exact four-body calculation 0. It has been shown 
that in 2D the second term is much more important than in 
3D. Without this contribution, the fluctuation contribution to 
the particle density, is divergent ll^ Is^ . The full 

derivative leads to a convergent particle density and hence an 
appropriate description of the BCS-BEC crossover. 
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III. STRONG COUPLING LIMIT: WEAKLY Here nis - 2m is the mass of the composite bosons. In the 

INTERACTING 2D BOSE CONDENSATE mean-field theory, the boson-boson coupling is a constant. 


While the mean-field theory predicts a simple and analyti¬ 
cal EOS in the entire BCS-BEC crossover, it does not capture 
correctly the interaction between the composite bosons in the 
strong coupling (BEC) limit. At yu = fic, the system undergoes 
a second-order quantum phase transition from the vacuum to 
the dilute BEC of bound states. In the grand canonical en¬ 
semble, the BEC limit corresponds to the regime - fXc + 0^, 
where the particle density n is vanishingly small. Alterna¬ 
tively, the chemical potential for composite bosons is given 
by 



(52) 


m 

This result is consistent with the previous calculation above 
the superfluid transition temperature j^ . However, it has 
been shown that for 2D bosons, the coupling is energy 
(chemical potential) dependent 114^1^ . that is. 


— - — In f 4 

.gB 47r , 


(53) 


Hb-2h + sb- 


(45) 


The BEC limit corresponds to //b —> 0 or, more explicitly, 
yUB/fiB ^ 0. 

Eirst, we show that the mean-field theory leads to a constant 
coupling between the composite bosons To this end, we 
derive the Gross-Pitaevskii free energy functional in the BEC 
limit Since the order parameter becomes vanishingly 

small for jj - fic + 0^, we can obtain a Ginzburg-Landau free 
energy functional of the order parameter held A(x), 


Ggl[A] = dx 






(46) 


The coefficients a, b, and c can be determined by the normal- 
state pair susceptibility Mo(iqi, q) which is given by (l33l) . Eor 
qh q^/(4OT) ^ fiB and jUb —> 0 ^, we have 


m ( \ 


Therefore, we obtain 


m »ZjUb 

a — b — -, c — -. 

47rsB 47reB 


(47) 


(48) 


The coefficient d contains the information of the interaction 
between the composite bosons. In the mean-field theory, it 
can be obtained by making the Taylor expansion of Qmf near 
A = 0. We obtain 


d = 


m 

4jtsI 


(49) 


Therefore, if we define a new condensate wave function 


ip(x) = 



(50) 


where ab is the boson-boson scattering length and y — 0.577... 
is the Euler constant. The constant coupling, ( |5^ . indicates 
that the BEC limit of the 2D mean-field theory corresponds to 
the Born approximation for four-body scattering in 2D jo^l . 
This is also true for 3D. However, in 3D, the Born approxi¬ 
mation already predicts a weak coupling .gB = ^nasImB with 
ab = 2a3d, and hence the 3D mean-field theory is qualitatively 
correct. 

Second, the incorrect boson-boson interaction can also be 
seen from the EOS. In the mean-field theory, the grand canon¬ 
ical EOS in the BEC limit can be expressed as 

^^mf(/^b) =(54) 
ion 

However, it is known that the grand canonical EOS of weakly 
interacting 2D Bose gases in the Bogoliubov theory is given 

byiH 


f^(/UB) = - 


mBfil 

Stt 


In 


_] 


1 

2 


(55) 


It was shown that the corrections beyond the Bogoliubov 
theory can be expanded in powers of the small parameter 
1/ln[4/(jUB'«BA^e^^^')]. The leading-order correction was 
presented in lIMIl . In the BEC limit /jb —> 0, the beyond- 
Bogoliubov corrections are vanishingly small in comparison 
to the Bogoliubov contribution. We expect that the Bogoli¬ 
ubov EOS, (l55l l. can be recovered in the BEC limit ^ibIsb ^0 
if we include the contribution from the Gaussian quantum 
fluctuations. If so, this allows us to determine the compos¬ 
ite boson scattering length by comparing our EOS with the 
Bogoliubov EOS, (l55l) . in the limit ^ibIsb —> 0. 

In the GPE theory, the grand canonical EOS is given by 


the Ginzgurg-Landau free energy reduces to the Gross- 
Pitaevskii free energy of a dilute Bose gas. 


^(Mb) 


mBfil 

8 n 


m + 






, ^B| |4 

v+—\v\ 


(51) 


where ^ - Hb /sb and the function f{0 is given by 


( 56 ) 
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f (0 = 




^OO ^(X 

I dx j 

Jo Jo 


dy In 


Here the dimensionless functions S, C, D, and are given by 

de r 

^U,y) = — 

Jo 27r Jo 


^ 2^2 y)C{x, y) + y'^Sjx, y)D{x, y) + TT^jx, y) ^ ^ C^x, y) + y'^D^jx, y) 

^ JlHx,y)+y^SHx,y) ^ JlHx,y) + y^!BHx,y) 


Six 

C(x, 

Dix, 

Tix. 


^ d0 r 

•"’'J. sj. 


dz 

) 

dz 




(£++^+)(£-+^_)l 

z+1 4 

(E, eJ 

(£+ H-£_)2 -H y2 


■^>‘1 SI 


4£+£_ (£+ + £_)2 + y 

i(±, > 

4U. 


1 1 
^ j (£+ + ^+)(£- + ^_) (£+ + £_)2 + y 2 
1 

,2 j_ „2 ’ 


, d0 r 

r^” d9 r°° 1/1 1 \ 1 

'’^^~Jo ^Jo '^^4(^^:^j(£+ + £_)2+y2’ 


4£+£_(£+ + f+)(£- + ^-) (£+ + £-)2 + y 

1/1 1 ' 


(57) 


(58) 


where the dimensionless variables x,y, and z are defined as x = q^/(4inBB), y = coIeb, and z - \JUniBB) and the notations 
and E± are given by 


= ;i(z + x + 2 V«COS0+1Ei = -y^(^)2T7. 


(59) 


We can show that /(^) is divergent at 4" = 0. To this end, we 
evaluate the functions S, C, D, and 'F at ^ = 0, which is 
denoted by the subscript 0. We have 

■^oix,y) = i In [(1 + x)^ + y^], 

1 y 

Su{x,y) = - arctan-, 

y 1 + X 

Jo L^(L^ + y2) (l 2 - 4xz)^^ 

r“ — 7x7 

Daix,y) = I dz - 

Jo L3(L2+y2)(L2-4xz)^^^ 

To{x,y) = r dz— - \ :■ (60) 

Jo (L2 + y2) xjlp- - 4xz 

Here we define L = z + \ + x for convenience. In the infrared 
limit, X ^ 0 and y —» 0, the above functions behave as 

JlQix,y) X, Soix,y)=^l, 

Coix,y)-^, £)o(x,y)^i, 

Toix,y)-]^. (61) 

For further analysis it is convenient to employ the polar co¬ 
ordinates X - p cos ip and y = p sin p. By making use of the 
Taylor expansion for the logarithm in (l57l i (see Appendix A), 
we find that at precisely ^ — 0, the function /(^) is divergent 
because of the infrared behavior JFq -h y^S^ ^ p^. We note 
that this kind of divergence does not exist in 3D. In 3D, the 


mean-field theory already predicts a weakly interacting Bose 
condensate in the strong coupling limit with a composite bo¬ 
son scattering length ob - 2 a 3 o a. The inclusion of the 
Gaussian contribution in the BEC limit leads to a modification 
of the composite boson scattering length from the mean-field 
value 2a3D to O.SSaso 

The divergence of the function /(^) at f —» 0 is not sur¬ 
prising. It is actually consistent with the Bogoliubov EOS 
(l55l l where the logarithmic term in the brackets diverges when 
Pb —> 0. Therefore, we expect that for ^ > 0, the function 

/(^) diverges as -In^ = ln(eB/pB)- To show this logarith¬ 
mic divergence, we separate the function /(^) into a divergent 
piece and a finite piece. The details are presented in Appendix 
A. The divergent piece is given by 

8 r'°° 

To capture the asymptotic behavior of this divergent piece for 
^ » 0, we find that it is sufficient to expand the denominator 
^ + y^S^ to the order D(f^) and approximate it as 

j?[2 + 3,2^2 ^ j(^x,y) ^ ml + y^Sl + Umo^x + (63) 


The explicit form of the function J?[i(x,y) is shown in Ap¬ 
pendix A. In the infrared limit p —> 0, we have m\ - 1. The 
neglected terms in the above ^ expansion lead to vanishing 
contributions for ^ ^ 0. Therefore, the infrared divergence in 
the limit ^ > 0 behaves as 


8 p/2 p 1/4 

d(p pdp - ^~ln—. 

71 Jo Jo p^+2^pcos<p + Pb 


( 64 ) 
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Thus we have shown that in the BEC limit 0, the Gaus¬ 
sian contribution Qqf behaves exactly like the logarithmic 
term in the Bogoliubov EOS, (BsT l. 

To obtain the composite boson scattering length a^, we 
need to determine the finite piece A, which can be defined as 

d = lim[/(0 + ln^]. (65) 


Using the definition of the fermion scattering length 020 , 


fiB 


4 


( 66 ) 


we obtain the composite boson scattering length 


flB - K:a2D, 



(67) 


A careful numerical analysis (see Appendix A) shows that A ^ 
-0.54. Therefore, we obtain 


K ^ 0.56. 


( 68 ) 


The energy density and pressure in the mean-field theory are 
given by 

Emf)'*) = ^^MF(/tMF) +/^MFH = EpG “ 2”®®’ 

7’mf(«) =-^^mfIAmf) = T’fG; (71) 

where Efg = nspl2 and Pfg = nepl2 are the energy den¬ 
sity and pressure of a noninteracting 2D Eermi gas with den¬ 
sity n, respectively. We see clearly from the pressure that the 
mean-field theory does not recover a weakly interacting Bose 
condensate in the strong attraction limit. 

To show that the chemical potential and the energy density 
suffer from the same problem, we define two dimensionless 
quantities 

^l + eB|2 E + nsB/2 

v=-, - - -. (72) 

£f r^FG 

In the mean-field theory, the solutions of v and R are in¬ 
dependent of the attraction strength in the entire BCS-BEC 
crossover; i.e.. 


This result is in good agreement with k ^ 0.56 from the 
exact four-body calculation and k ^ 0.55(4) from the 
EOS predicted by the diffusion Monte Carlo simulation 1^ . 
We also notice that the pole approximation of the Gaussian 
quantum fluctuations with a dimensional regularization of the 
UV divergence in the BEC limit predicted an analytical result 
4 = -1/2 and hence k = 1 l(2^l^e^l^) ^ 0.55 III. 


IV. BCS-BEC CROSSOVER 


In this section, we study numerically the EOS in the en¬ 
tire BCS-BEC crossover. The determination of the grand 
canonical EOS is simple. The grand potential Q(ju) can be 
obtained by performing the numerical integration in ( l38l l for 
—ebI' 2 - < fi < -Hoo. The BEC and BCS limits corresponds 
to /i —> —bbI2 and ji —> -too, respectively. In this work, we 
are interested in the canonical EOS for a homogeneous 2D 
Eermi gas with fixed density n. This enables us to compare 
our results with recent quantum Monte Carlo calculations of 


the energy density 


of the local pressure 1137 


4711 and experimental measurements 
m- 'or convenience, we define the 


Eermi momentum ^f and the Eermi energy ep for a noninter¬ 
acting 2D Eermi gas with the same density n. They are given 


by ^F = yflmi and sf = nnlm. The BCS-BEC crossover is 
controlled by the dimensionless ratio a = Sb/sf or the gas 
parameter 


7] - \n{kpa2D)- (69) 

The BCS and BEC limits correspond to rj ^ -i-oo and 77 ^ 
- 00 , respectively. 

In the mean-field approximation, we have n ^ tiyipiM), 
which gives rise to the mean-field results of the chemical po¬ 
tential and the pairing gap lHUlU], 

£b I 

RmfOi) - sf —Amf(«) = ■\J2eBEF- (70) 


vmf - 1, Rmf - 1- (73) 


On the other hand, the Bogoliubov theory predicts that the 
canonical EOS of a 2D Bose gas is given by ll76l - [^ 


47rnB 1 



E - -ubEb + 


2nn\ 1 



, (74) 


where hb - n 12 is the density of tightly bound bosons. There¬ 
fore, we expect that in the BEC limit (77 —> - 00 ) the solutions 
of V and R behave asymptotically as 

''~2in(^)- 277’ ^~2in(^)_2;7’ 

where k ^ 0.56 from the exact four-body calculation or 
from our study in Sec. HII] These results indicate that v and R 
become vanishingly small in the BEC limit. We note that the 
use of the Bogoliubov EOS (l74l) requires that the parameter 
1/ln[l/(? 7 Ba|)] is sufficiently small or 77 —» - 00 . The correc¬ 
tions beyond the Bogoliubov theory was studied in Refs. izi- 
1^ . On the other hand, in the BCS limit (77 —» -Foo), the 
pairing gap A becomes vanishingly small and hence the GPE 
theory becomes equivalent to the particle-particle ladder re¬ 
summation lfl4l - [T^ . Therefore, in the BCS limit, the GPE 
theory naturally recovers the perturbative EOS of a weakly 
interacting 2D Eermi gas up to the order 0 ( 1 / 77 ^). The per¬ 
turbative EOS of a weakly interacting 2D Eermi gas is given 

byimm 


V 

R 


, 1 r-Fl- 21 n 2 

1 - + - - - 

77 77 ^ 

1 7 H- 3 / 4 - 21 n 2 

1 “ H" 9 

?7 


+ 0 


+ GI-I. 


(76) 


Therefore, we expect that v and R approach unity asymptoti¬ 
cally for rj +oa. 
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In the GPF theory, the chemical potential ^ is determined 
by solving the full number equation 

n = «mf(ju) + «gf(ju). (77) 

Then we can determine the energy density E{n) - Qmf(A() + 
Ggf(/^) + pw and the pressure P{n) - - ^Gpip)- The 

Gaussian contribution ncpifi) can be worked out analytically 
but it is rather tedious. In practice, we start from the grand po¬ 
tential Q.(jj.) - Qmf(A') + To determine the chemical 

potential //, we calculate the energy density as a function of yu; 
i.e., E{fi) - Q(yu) -I- fin. We search for the maximum of E(ju), 


which gives rise to the solution of the chemical potential for 
the given density n. Meanwhile, the energy density and the 
pressure for the given density n are determined. To perform 
the numerical calculation, it is convenient to use the dimen¬ 
sionless variable v. The mean-field contribution to the grand 
potential is Q.mf(m) - -v^Ffg- The Gaussian contribution to 
the grand potential can be expressed as 

^gip) = .g(i')£FG, (78) 

where the function g(v) is given by 


J 


.g(v) = 



2 2 t) + fS{s, t)D{s, t) + 2T^{s, t) 

^ J^{s, t) + f) 


+ 16a^v^ 


C\s, i) + fD\s, t) 
t) + fi&is, t) 


The dimensionless functions S, C, 29, and T are now defined as 


^(s, t) 
S(s, t) 
C{s,t) 
29(s, t) 
r{s,t) 





(£++^f)(£-+^_)1 

2u + a 4 

[E, eJ 

(£+ -F £_)2 -F f2 


(£+ 


do p 

Jo 27 r Jo 

Jo 27rJo £_ 


1 (£++^+)(£.+f) 


-F £_)2 -F f2 
1 


1 


(£++^+)(£-+^_)(£+ + £_)2 + f2 


de r°° i i 

Jo 2^ Jo 4£+£-(£+ + ^+)(£- + ^-) (£+ + E-)^ + fl ’ 


lo 2nJo 4Uf £-/(£+ + £-)2 + f2’ 


(79) 


(80) 


where the variables s - q^/iSmsp), t = and u - 

k^/(2msY^). Here the notations and E± are given by 

- u + s ±2yuscos6 - V + 

E^ = ^(^j 2 + 2 av. (81) 

Using the function g(v) we have defined, we can express 
the dimensionless quantity E as 

E(fl) -F hlSB ^ 

Riv) = - ^ - = -v" + .g(v) + 2 v. (82) 

^FG 

The physical results of v and R correspond to the maximum 
point of the the function R{v) in the range 0 < v < 1. In the 
mean-field theory, we neglect the Gaussian contribution g(v) 
and hence R(v) ^ -v- + 2v. The maximum of the function 
R{v) gives the results v = 1 and R - 1, which are precisely the 
mean-field predictions, (17^ . Including the Gaussian contri¬ 
bution g(v), the maximum of R(v) will be modified since the 
function g(v) depends explicitly on the interaction strength a 
or the gas parameter 77 . In Fig. [1] we show the curves of the 
function R(v) for several values of the gas parameter 77 . We 


find that the quantum fluctuations become more and more im¬ 
portant when the attraction strength increases. 



FIG. 1: (Color online) Curves of the function R(v) for various values 
of the gas parameter 77 = ln(feFa2D)- For comparison, we show the 
mean-field prediction = -v^ + 2v by the dashed line. 
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FIG. 2; Evolution of the chemical potential yu in the BCS-BEC 
crossover. We show the quantity v = (ju + e^I2)Ie^ as a function 
of the gas parameter rj = \n{k^a 2 o)- The mean-field prediction is 
represented by the dashed line. 



FIG. 3: The order parameter or pairing gap A (divided by ep) as 
a function of the gas parameter 77 = ln(^:pa 2 D)- The dashed line is 
the mean-field prediction. The dashed-dotted line corresponds to the 
prediction with GMB effect in the weak coupling regime (77 > 2). 


In Fig. [2] we show the evolution of the chemical potential ji 
or, explicitly, the quantity v - {^ + £B/2)/eF in the BCS-BEC 
crossover. We find that v ^ 1 in the BCS limit and v ^ 0 
in the BEC limit, in agreement with our general expectation. 
The order parameter A is shown in Eig. [3] We find that the 
inclusion of the quantum fluctuations leads to a large suppres¬ 
sion of the order parameter in the strong coupling regime. At 
weak coupling, it was shown that the induced interaction or 
the Gor’kov-Melik-Barkhudarov (GMB) effect 19^ leads to 
a suppression of the critical temperature and hence the pairing 
gap A by a factor of 1/e In Eig. [3l we also show the 
prediction with the GMB effect in the weak coupling regime 
(77 > 2). Obviously, the current GPE theory does not take into 
account the GMB effect. 

In Eig. m we show the evolution of the energy density E 



FIG. 4: (Color online) Evolution of the energy density E in the BCS- 
BEC crossover. The quantity R - {E + 77£b/2)/£'pg is shown as a 
function of the gas parameter 77 = ln(fepa 2 D)- The dashed line repre¬ 
sents the mean-field prediction. The (blue) circles and (red) squares 
represent the predictions from the diffusion Monte Carlo simula¬ 
tion 0 ] and the auxiliary-field Monte Carlo simulation respec¬ 
tively. The bottom-left dashed (green) line represents the Bogoliubov 
EOS of a weakly interacting 2D Bose gas with the boson scattering 
length Ob = 0.56a2D [see Eq. (Hi)]. The top-right dashed (purple) 
line shows the EOS of a weakly interacting 2D Eermi gas [see Eq. 

nil]. 



FIG. 5: (Color online) Evolution of the pressure F in the BCS-BEC 
crossover. P/Ffg as a function of the gas parameter 77 = ln(kFa 2 D) is 
shown. The mean-field prediction is represented by the dashed line. 
The (blue) circles with error bars are the experimental data taken 
from 


or, explicitly, the quantity R - (E + 77eB/2)/£’FG in the BCS- 
BEC crossover. In the BEC limit (77 —> - 00 ), our result ap¬ 
proaches the Bogoliubov EOS, (|74|, of weakly interacting 2D 
Bose gases with the boson scattering length ab - 0.56a2D- In 
the BCS limit (77 —> - 1 - 00 ), our result tends to the perturbative 
EOS, (|76] |, of weakly interacting 2D Eermi gases. The energy 
density was com pute d recently by using the diffusion Monte 
Carlo simulation and the auxiliary-field Monte Carlo sim- 
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FIG. 6: (Color online) Evolution of the contact C in the BCS-BEC 
crossover. This figure shows (C - Cmf)/^f as a function of the gas 
parameter r] = ln(/:Ffl 2 D): where Cmf = a/2 is the mean-field predic¬ 
tion. The dashed (blue) line and the dash-dotted (red) line represent 
the predictions from the diffusion Monte Carlo simulation 14611 and 
the auxiliary-field Monte Carlo simulation El , respectively. 


ulation El. In Fig. m we also show these Monte Carlo results 
for comparison. Even though our theory recovers the correct 
BCS and BEC limits, there exists a slight deviation between 
our theoretical prediction and the Monte Carlo results. This is 
not surprising because the GPE theory, which considers only 
the Gaussian pair fluctuations, is not an exact treatment. Some 
many-body effects we have not taken into account in the GPE 
theory may account for this disagreement. Eirst, the current 
GPE theory does not consider the GMB effect 1^ . which 
leads to a suppression of the pairing gap A by a factor of 
1/e at weak coupling El- The inclusion of this effect may 
lead to a slight suppression of the energy density and a faster 
convergence to the EOS, (fThl) . of weakly interacting 2D Eermi 
gases. The GMB effect may also be important in the crossover 
regime (roughly -0.5 < tj < 2). Second, in the GPE theory, 
we consider only the Gaussian pair fluctuations. The contri¬ 
butions from the non-Gaussian quantum fluctuations (beyond 
quadratic order in 0 and (p*) may be important to make for a 
better agreement with the Monte Carlo results in the crossover 
regime. 

As we have mentioned, the most important thermodynamic 
quantity which shows the significance of the quantum fluc¬ 
tuations is the pressure P. The mean-field theory predicts 
P = Pgf for arbitrary attraction strength. In the GPE theory, 
we have 


P(n) 

Pfg 


2 v-R, 


(83) 


where v and R have been determined by searching for the max¬ 
imum of the function R{v). Therefore, the pressure depends 
explicitly on the interaction strength. In Eig. |5] we show 
the evolution of the pressure or, explicitly, the ratio PI Pfg in 
the BCS-BEC crossover. Recent experiments on the quasi-2D 
Eermi gases across a Eeshbach resonance have measured the 


local pressure at the center of the atom trap at sufficiently low 
temperatures EIEIj which can be regarded as the ground- 
state pressure of a homogeneous 2D Eermi gas in the BCS- 
BEC crossover. In Eig. |5l we also show the experimental data 
reported in El . Except for the deep BCS regime (77 > 3), 
our theoretical prediction is in good agreement with the ex¬ 
perimental measurement. The observed high pressure in the 
deep BCS regime could be attributed to the mesoscopic na¬ 
ture of the experimental system: In the deep BCS regime, the 
scattering length 1220 becomes larger than the cloud size and 
hence the interaction is effectively suppressed El . On the 
other hand, it has been argued that the temperature effect may 
also be crucial to understand the observed high pressure in the 
deep BCS regime El . In the future, it is necessary to study 
the finite-temperature effect in the current GPE theory. 

Having determined the EOS, we can calculate the contact 
C, which is a powerful quantity to relate the ener gy, p ressure, 
and the microscopic momentum distribution 19511 . In 2D, 
the contact C can be defined as 

C _ 1 d{E/E,G) 

id A J ■ VO'+J 

4 dri 

After some simple manipulation, we obtain 


C /J. 

E 

1 ^ 

E ] 

fiF 

Efg 2 

[Pfg 

EfgI 


Using the mean-field result Cmf/^p = a 12, we can show that 


C - Cmf 

kl 


■R. 


( 86 ) 


In Eig. | 6 ] we show the quantity (C - Cmf)//^p in the BCS- 
BEC crossover. We find that this difference is quite small in 
the entire BCS-BEC crossover and is peaked around 77 ^ 0.7, 
which agrees with recent quantum Monte Carlo results ES 

0 . 


V. SUMMARY 


The lack of a weakly interacting Bose condensate in the 
strong attraction limit is a longstanding problem for the the¬ 
ory of BCS-BEC crossover in two-dimensional Eermi gases. 
Especially, the mean-field prediction for the pressure in the 
BCS-BEC crossover shows the inadequacy of the mean-field 
theory in 2D. The inadequacy of the 2D mean-field theory can 
be understood from the fact that the Born approximation for 
four-body scattering in 2D predicts an incorrect form of the 
composite boson coupling. In this work, we showed that this 
problem can be solved by including the contributions from the 
Gaussian quantum fluctuations. In the BEC limit, the missing 
logarithmic dependence on the boson chemical potential and 
hence the boson-boson interaction is naturally recovered by 
the quantum fluctuations. We determined the composite bo¬ 
son scattering length as - 0.56a2D, in good agreement with 
the exact four-body calculation and recent quantum Monte 







12 


Carlo results. We calculated the chemical potential, the en¬ 
ergy density, the pressure, and the contact for a homogeneous 
2D Fermi gas in the BCS-BEC crossover. Our theoretical pre¬ 
dictions are in good agreements with recent quantum Monte 
Carlo results and experimental measurements. 

In the future, it is necessary to consider more many-body 
effects to explain the slight discrepancy between our theoret¬ 
ical prediction and the quantum Monte Carlo results, such as 
the GMB effect and the non-Gaussian fluctuations. In the 
BEC limit, an exact low-density expansion for the compos¬ 
ite bosons iH could also exist in 2D. It is also interest¬ 
ing to extend the present theoretical approach to the finite- 
temperature case and the spin-imbalanced case. The inclu¬ 
sion of the Gaussian fluctuations may provide better predic¬ 
tions for the Berezinskii-Kosterlitz-Thouless transition in the 


2D BCS-BEC crossover an d the phase structure of spin- 
imbalanced 2D Eermi gases 11^ . 
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Appendix A: Counting the infrared divergence of the function /(O 


Using the Taylor expansion ln(l - a) = - a" jn, we can express the function /(^) as 


m 


4 r“ r“ J[C + y^B!D + 2T^ r“ r 

= Oo ‘'Go - -1 ‘'■‘Jo 

O ^ >'2/2—2 ^ / \ j^OO /^Cx 


(^C + /SD -H + /2)2)"-* 


To analyze the infrared divergence for f ^ 0, we expand the quantities and S in the denominators in powers of 

Jl(x,y) = JlQ{x,y) + V ^Jl„(x,y), S(x,y) = So{x,y) + V ^S„(x,y), 

^ nl ^ nl 


(Al) 


(A2) 


n=l 


«=1 


where and = d''S j The expansion coefficients Ji,, and can be evaluated to arbitrary order by 

using Mathematic a. Here we list the results for and ‘B\. We have 


= 

Si = 


i ‘ 

-r 


1 


L2 + y2 

2L 


1 H- 


%xz 


dz 


(L2 - 4xz) 
1 


3/2 


£2 + ^2 ^^2 _ 4 ^^ f /2 


f 


f 


dz- 


dz- 


2 L 2 


(L 2 + y2^ \ yJjJ- - ‘Xxz 


-1 , 


2L 


1 - 


(L2 -i-y2)2 ( yjjj. _ 


(A3) 


where L = z+\ H-xas defined in the text. In the infrared limit x,y 0, we have iTli ^ 1 and Si —> -1. The expansion of the 
quantity + y'^& takes the form 


OO n 


Jl^(x, y) + y^S^ix, y) = ^(x, y) H- y^S^ix, y) + ^ T ^ ,,, [^k(x, y)Jln-k{x, y) -H y^Sk(x, y)S„-k{x, y)]. (A4) 


n=\ k=0 


k\(n-k)\ 


Eor further analysis, it is convenient to use the polar coordinates x - p cos ip and y = p sin ip. At exactly ^ = 0, we have 
oZIq + y^S\ ^ in the infrared limit p —> 0. To capture the leading asymptotic behavior, we find that it is sufficient to 
approximate tfie quantity -i- y^S^ as 


J^{x,y) +y^S^{x,y) ^ J(x,y) = Jll(x,y) + y^Sl(x,y) + 2^Jlo{x,y)Jli{x,y) + ^^Jl1(x,y). 
In the infrared limit p ^ 0, the function JTix, y) behaves as 

J{x, y)^p^ + 2(p cos (p + ^^. 

The other contributions we neglected in approximation (IA5b behave in the infrared limit as 

O \ /CO \ /CO 


p^+ p+ 


\/2=l / 


V/i=2 


Vn=3 


(A5) 


(A 6 ) 


(A7) 
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The above terms lead to vanishing contributions in the limit ^ ^ 0. One can prove this observation by carefully analyzing the 
infrared behavior of the following integral 


r p"' 

tmn — I pdp 2 _i_ or ~i~T2\n' 

Jo + 2^pCOS(/3 + 

The properties of the integral I„,„ can be summarized as follows: For m > 2{n - 1), the integral is finite; for m - 2{n - 1), it 
diverges as I,„„ -In^; and for m < 2{n - 1 ), we have ~ jy'^ 2 n- 2 -m ^ ^ q 



FIG. 7: The quantity f(() + ln(0 as a function of (in the range 10^<^<10^. In the calculation we use Eq. (A 16) for f((). 


In the infrared limit p —> 0, the second term in the expansion, (lAll) . behaves as 
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It vanishes in the limit ^ ^ 0. The third term in the expansion, (lAll i. behaves as 
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A careful analysis shows that this term leads to a finite contribution in the limit ^ > 0. The nonvanishing contribution from the 
k - n terms can be expressed as 
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Next we analyze the first term in (lAlb . which develops the logarithmic divergence. In the infrared limit, it behaves as 
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Therefore, we can separate the above contribution into two pieces. The finite piece is given by 
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The divergent piece is given by 


In the infrared limit, this piece behaves as 
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(A15) 


Therefore, it exactly develops the asymptotic behavior /(^) ~ - In f + d for ^ ^ 0. Summarizing the nonvanishing pieces 
(All), (A13), and (A14), we find that to capture the logarithmic divergence and determine the finite term A, it is sufficient to 
approximate the function /(^) as 
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(A16) 


In Fig. |2]we show the numerical result of f(^) + In(^) in the range 10 ® f < 10 It is clear that in the limit ^ ^ 0, it converges 
to a constant. Thus we determine A ^ -0.54. ._ 
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